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We present the energy density and flux distribution of a heavy quark moving through high temperature QCD 
plasmas and compare them with those in the strongly coupled Af = 4 SYM plasma. 

The Boltzmann equation is reformulated as a Fokker-Planck equation at the leading log approximation and 
is solved numerically with nontrivial boundary conditions in momentum space. We use kinetic theory and 
perform a Fourier transform to calculate the energy and momentum density in position space. The angular 
distributions exhibit the transition to the ideal hydrodynamics and are analyzed with the first and second order 
hydrodynamic source. 

The AdS/CFT correspondence allows the same calculation at strong coupling. Compared to the kinetic theory 
results, the energy-momentum tensor is better described by hydrodynamics even after accounting for the differ- 
ences in the shear viscosities. We argue that the difference between the Boltzmann equation and the AdS/CFT 
correspondence comes from the second order hydrodynamic coefficient which is generically large compared 
to the shear length in a theory based on the Boltzmann equation. 



1. Introduction 



When a heavy quark moves supersonically through plasmas, the energy and momentum are redistributed, 
producing a Mach cone and diffusion wake. Originally, the study of the Mach cone structure was motivated by 
an attempt to understand the two particle correlation functions measured in the heavy ion collisions. Although 
the Mach cone failed to explain the experimental observation, it is a simple way to study the medium response 
to an energetic particle. 

Using the AdS/CFT correspondence, the stress-energy tensor in the strongly coupled SYM plasma due to 
a heavy quark have been calculated pfl In this work, we study the Mach cone produced in the high 
temperature QCD plasmas at weak coupling and compare it with the strong coupling results given by the 
AdS/CFT correspondence. 



2. Linearized Boltzmann Equation at Leading Log 



We start with the Boltzmann equation 



where the particle distribution function is linearized around the equilibrium, /(t,x, p) = Up + Sf{t,x,p) with 
Up — l/(eP/"^ — 1). The collision term is given by 



we consider only t-channel 2 — > 2 scattering. Two hard particles with momenta of order T are scattered into 
two others, exchanging soft particles with momenta of order gT, where g is the weak coupling constant. The 
particles are scattered with angles of order g. With the definition 6f{p) = np{l + np)x{p), the collision integral 



(9t + Vp.ax)/(i,x,p) = C[/,p], 



(1) 




2 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



can be simplified and expressed as a variational form ^ [21 S] 



d f d 
{dt + Vp ■ d^)6f{t,x,p) = TfiA-Q-i \np{l + np)—^ 



Sf 



np{l + Up) 



gain terms , 



Svr 



T 

mo 



Without gain terms, the equation can be reformulated as a Fokker-Planck equation 

{dt + Vp • d^)df{t, X, p) = -fiA (1 + 2np) P • + r/i^V2 Sf 



(3) 



(4) 



and the motion of the particle excess, Sf, can be described by the Langevin equation with the drag and random 
kicks 



dp' 
~dt 



-fiA{l + 2np)f + f (t) , {C{t)e{t')) = 2TfiAS''S{t - t') . 



(5) 



In this diffusion process, the excess particles lose their energy and momentum to the plasma. However, there 
are work and momentum transfer on the particles per time, per degree of freedom, per volume 



dE 
dT 
dP 



= -T^A 



^Ml + n,)p.- 



(27r)3 ' "'dp- 

In terms of these energy and momentum transfer, the gain terms are given by 



gain terms — 



)_d_ 
dp 



p^np{l + Up) 



dE 
dtT 



1 



d_ 
dp 



np{l + Up) 



dP 



(6) 



(7) 



where we defined = / ^.,^3 ^p(l + Tip) = Using the Boltzmann equation with the gain terms, it is easy 
to show that the lost energy and momentum of the excess particles are exactly compensated by the gain terms. 
Therefore, the total energy and momentum are conserved during the evolution. 

In order to get the boundary condition, we consider the excess of soft gluons within a small ball of radius 
Ap ^ gT centered at p = 



P=^P d^p 



p=0 



^np{l + np)x{p) 



x{0)Ap 



(8) 



which is vanishing at the weak coupling limit. Since it is easy to emit or absorb a soft gluon, we choose the 
absorptive boundary condition 5j 



X(P)I 



p-!-0 



= 0. 



(9) 



As a result, there is flux (dx/dp) at p = and the particle number changes. 

Using the linearized Boltzmann equation with the absorptive boundary condition, we determined the transport 
coefficients including the shear viscosity rj and the second order hydrodynamic coefficient [Ij 



r? T 
— = 0.4613 — 
sT fiA 



= 6.32- 



sT 



(Boltzmann) . 



(10) 



which will be used in the Mach cone analysis in the following. In comparison, we provide the corresponding 
coefficients for the AdS / GET correspondence j6l [7] 



rj 1 
'sf ~ AttT 



2 -In 2 
2ttT 



(AdS/CFT) 



(11) 



^The dimension of the adjoint representation is dj^ = — 1. The Casimirs of the adjoint and fundamental representations are 



Ca = Nc and Cp = (N^ - l)/2Nc. The Debye mass is m^, = 2g'^CA f 
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3. Mach Cone 

We simulate a Mach cone by solving the linearized Boltzmann equation with a source term 

{dt + Vp • ax) 5f{t, X, p) = C[f, p] + S{t, X, p) . 



(12) 



In the presence of a heavy quark moving with a constant velocity v, the particles in equilibrium are scattered, 
producing the source around the quark, 5(t,x,p) = S'(p)(5'^(x — -vt). The source term is given by another 
collision integral 



5(p)--^ / \M\''{2Tiy5^{P+K-P'-K') 

2 Jk,p',k' 



/(p)[l + /(p')]<53(k-k^)-/(p')[l + /(p)]<5-^(k'-kH) 



, (13) 



where we used (k) — (27r)^(5'^(k — kH)(5'^(x — vi) as the distribution function of the heavy quark. Especially, 
when the quark moves with the speed of light, w = 1, it is simplified at the leading log order 



l + 2np\ (l + 27ip) 



T 



T 



p V 



g'^Cpnn? 
8^ 



^log 



(14) 



Our strategy to get a Mach cone is the following: (i) perform a Fourier transform of the Boltzmann equation 
and solve it in Fourier space, (ii) calculate the stress-energy tensor in Fourier space using kinetic theory, (lit) 
perform a Fourier transform of the tensor back to position space to get the distributions. We assume the 
heavy quark moves along z in the lab coordinate system (x, y, z) and introduce the Fourier coordinate system 
(x',?/',z'), where z' points along the Fourier momentum k. For simplicity, we use the rotational symmetry 
around z and let k lie on zx plane. Then the source in the Fourier coordinate system becomes 



5(p) = V47r 



_JJrF_ 
2dA^l 



-rtp(l + Up) 



1 + 2nr 



T 



(1 + 2np) 



cos6 Hifi{p) + sin9 _ffi,i(p) 



,(15) 



where 9 is the angle between z and k and Hi^„i{p) is the real spherical harmonics defined in Fourier space [4] 



ffo,o(p) 



47r 



The Boltzmann equation in Fourier space is given by 

{-iu + iwp ■ k)5f{u, k, p) = C[Sf, p] + 27rS'(p)(5(a; - v • k) , 



(16) 



(17) 



We solve it for (5/(cj,k,p) = 2Tr6{uj — v • k)?ip(l + np)x(i^,k,p) and calculate the stress-energy tensor using 
kinetic theory 



ST"^{uj,k)^2dA I p^<5/(c^,k,p) 



(18) 



By Fourier-transforming back to position space, we get the energy and momentum density distribution 

3-''"*+^'"^<5T™(w,k), 



5r™(t,x) = 2d A f 

6T°'it,x) = 2d A 
dT"^t,x) = 2d A 



-ic<ji+«k-x 



sin e ST°''' (w, k) + cos 6 ST°''' (w, k) 
cos ST^'-^' (w, k) - sin 6 ST"^' {lo, k) 



(19) 



where the integrands inside the square brackets are the real momentum ST'^^{uj,'k) and (5r°^(a;,k) in terms of 
Fourier components. The numerical results of the energy and the momentum density distribution are shown 
in Figure [l] In comparison, we present the AdS / CFT results on the same scale in Figure [2] 
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-R[A£(R)] 




5 



Figure 1: (a) The energy density distribution, _R [A£'(R)]/[/i_F(CsSr)/(477/3)] and (b) the energy flux distribution, 
R [A5'(R)]/[^F(CsSr)/(4r;/3)] given by the Boltzmann equation. They show a smooth transition from the free streaming 
to the hydrodynamic regime. We can see a Mach cone and a diffusion wake (in the flux distribution) behind the heavy 
quark. (R = a;TX + s^z) 



fl[A£(R)] 




5 



Figure 2: (a) The energy density distribution, i?[A£;(R)]/[7r^r^^A/(l - v'^)/2] and (b) the energy flux distribution, 
ii[AS'(R)]/[7r^r^'\/A/(l — «^)/2] given by the AdS/CFT correspondence. They are singular near the heavy quark at 
origin. We can see a Mach cone and a diffusion wake (in the flux distribution) behind the heavy quark. (R = sitx + xlz) 



In order to analyze the structure, we define the angular distribution with an angle Or measured relative to z 
on zx plane. At a fixed distance R = Ixxx + Xizj from the heavy quark, the energy and the momentum density 
are given by 

27ri?2 sin6l;^ Ai;(R) , 

27ri?2sin6lfl R- AS(R) , (20) 

where l^E = and A^i = (5T°\ The results are shown in Figure [s] where distances (i?) are measured 

compared to the shear length, r]/{sT). Generally, the Boltzmann equation and the AdS/CFT correspondence 
give similar distributions with a Mach cone and a diffusion wake, except for 71 = 1. At short distances, the 
AdS/CFT solutions are dominated by the zero temperature contribution which is absent in the Boltzmann 
formulation. 



4. Hydrodynamics 

Hydrodynamics is an effective theory at long distances and is given by the conservation of the energy- 
momentum tensor 

Trydro = HT) + V{T)]u^u^ + V(T)g'^^ + i:^'' , (21) 



dER, 

d^R 
dSR 
d^R 
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Figure 3: The angular distribution of (a) the energy density [dER/d^Rj/^p and (b) the energy flux [dSR/d^Rj/^xp at 
distances TZ = 1, 5, 10, 20, and 40. The Boltzmann results are compared with the AdS/CFT results. The distance from 
the heavy quark is measured compared to the shear length, rj/{sT). At 7?. = 1, the AdS/CFT results are dominated 
by the zero temperature contribution which is absent in the Boltzmann formulation, hf for the Boltzmann equation 
corresponds to TTT^y''X/{l — v^)/2 for the AdS/CFT correspondence. 

where the first two terms on the right hand side correspond to the ideal hydrodynamics and the last term is 
the dissipative part. The heavy quark perturbs it away from equilibrium 



e{t, x) ~ Co + Se{t, x) , u^{t, x) ~ (1, u{t, x)) . 



(22) 



where 5e and u are the first order disturbance. In the linearized conformal limit, the dissipative part is given 
by the first and second order derivative of the fields [7] 



tt'"' = -277(9''it'') - 2rjT^{d''d''\nT) 



(Static) , 



(23) 



where ( • ■ ■ ) denotes the symmetric and traceless spatial component. The second order term can be replaced by 
the leading order relation, making it as a dynamical equation 



(Dynamic) 



(24) 



With the constituent relation, we solve the equation of motion in presence of the external force (minus the drag 
force) 



7mcro ' 



In this case, the force -Fmicro = dp^/dt is given by 

K^^cro = / P'Sit, X, p) = ^1FV''5''{^ - vt) 

J D 



(25) 



(26) 



where v'^ — (w^, v). 

At long distances from the quark, the stress-energy tensor is described by Tl^y^^^ and at short distances, we 
have to consider correction r'^" which is localized near the quark 



hydro 

Then the equation of motion in Fourier space becomes 



(27) 



(28) 
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Table I: The correction to the source, [i?] — \j3~^ 





01 (0,0) [R] 




<^2(0,0) [R^] 


Bohzmann 








0.48 


AdS/CFT 


-1 


0.34 


-0.33 



and the correction part acts as a source. In the second order, r'-' can be written with three symmetric tensors 
consisting of v and k 



2tt^f 5{uj — V • k) 



y^yJ _ l„2^y ) 0i(w,fc^) 



(29) 



where (l)i{uj, k"^) (i — 1, 2, 3) is analytic near a; = k = 0. Since r*-' is localized, we can expand it for small w and 
k 



ifcV*^ = 2TT^iF <5(w - V • k) 



3 lau) 



idbJ 



+ Oik^). (30) 



We determine r^"^ by comparing the numerical solution T^'' with the constituent relation T^\^^. By comparing 
order by order, all of the coefficients are numerically determined. Table |l] shows the result where the coefficients 
are measured compared to the shear length unit [i?] = [(4?7/3)/(CgST)]. Note that 4>i{0, 0) and ^"^-g"'"^ vanish for 
the Boltzmann equation. Since the source S'(p) in Eq. (fTsl) does not have i?i,±2 components which correspond 
to the tensor \y^v^ — , there is no (/)i(w, k'^) in tlieBoltzmann case. 

With the corrected force F^^^^^ — ikW^^ = 2T:fipS{uj — v • k)[v (/)„(w, A;^) + ik (/)fc(w, fc^)] [2], we solve the 
equation of motion with v — 1 



- iuj 6T°° + ik 6T'^ ^ 27r/iF '5(w - V • k) , 



ja;(5r°^' +«c^(fc)fc5T"" + r^(a;)fc^5T"^' = [cos 6i + i/c 0^] 27r^F ^(^ - v • k) 



—iu! ST' 



Ox' 



D{uj)k^ ST^""' = sm9(j)y2n^F6{uj -V -k) , 



(31) 



where <f>^(ui,k'^) and (/)k{uj,k'^) are expanded around a; = k = depending on the order and c^(fc) = = 1/3 
is the speed of sound, Ts{uj) = = (4?7/3)/(sr) is the sound attenuation length, and D{uj) = D = ri/{sT) is 
the diffusion constant in the first order hydrodynamics. In the second order static hydrodynamics (Eq . (|23[ )), 
we have same constants except c^(fc) = c^(l + TTrFsfc^) while the exception for the dynamic case (Eq. (|24[)) is 
rs(a;) = rs/(l — ir^^uj) and D(uj) ~ D/{i — iTT^uj). The hydrodynamic solutions are given by 



,00/, , _ ^ + fc cos 6> - fc^ [Fg (a;) + 0^] 
^2 _ c2(fc) fc2 iV s{uj) ujk^ 



„, i \uj cos (l)y + c^{k)k] ~ kuj<pk. 

61 [0J,k) — 

(5r°^'(w,k) 



2'KiJLp5{(jj — kcos9) , 



^2 _ c2(fc) fc2 jFs(a;) ajfc2 
i sin 9 d 



2'K^p5{uj — k cos ( 



w + iD{uj) k^ 



2'K^p5{uj — fc cos ( 



(32) 



In Figure |4j [5j we compare hydrodynamic results to the full theory. The difference between the static and the 
dynamic solution is shown in Figure [6j 



5. Discussion 



We presented the energy density and fiux distribution given by kinetic theory based on the Boltzmann equation 
at weak coupling, and compared them with those given by the AdS/CFT correspondence at strong coupling. 
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Figure 4: The angular distribution of (a) the energy density [dER/d^R]//iF and (b) the energy flux [dSR/d^R]//iF given 
by the Bohzmann equation at distances TZ — 10, 20, and 40. The Bohzmann results are compared with the first order 
and the second order static hydrodynamics. 




Figure 5: The angular distribution of (a) the energy density [dEH/d6'R]/[7r7\/AT^/2] and (b) the energy flux 
[dSR/d6lR]/[7r7VATV2] given by the AdS/CFT correspondence at distances 7^ = 5, 10, and20. The AdS/CFT resuhs 
are compared with the first order and the second order static hydrodynamics. (7 — 1/Vl — v^) 



Both of them show a Mach cone and diffusion wake at long distances but their beliaviors are distinguishable by 
the approach to the hydrodynamic limit. 

In tfie Boltzmann equation, the second order static hydrodynamics (Eq. ( |23[ )) works better than first order 
for TZ > 40, but it is unnatural for TZ < 10 due to the high frequency behavior. In principle, the second order 
dynamic solution (Eq. (24)) is better than static one at long distances (see Figure [6] (a)). However, the dynamic 
solution has higher order terms and develops spurious shocks for TZ < 20. 

In the AdS/CFT correspondence, even the first order hydrodynamics is good for TZ > 5. Only the sound 
wave has corrections of the second order static hydrodynamics and is dominant in the energy distribution. At 
TZ = 5 in the energy density, we can see the improvement given by the second order theory. Both the static and 
the dynamic solution work very well in Figure [6] (b) . 

In comparison with the AdS / CFT correspondence, the Boltzmann equation produces a smooth transition from 
the free streaming to the hydrodynamic regime. The stress-energy tensor given by the AdS/CFT correspondence 
is significantly better described by hydrodynamics in the sub-asymptotic regime. Although they have different 
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Figure 6: (a) The angular distribution of the energy density, [dER/d^Rj/z^p given by the Boltzmann equation at distances 
TZ = 20, 40. The Boltzmann results are compared with the second order static and dynamic hydrodynamics. The dynamic 
solution develops spurious shocks for TZ. < 20. (b) The angular distribution of the energy density, [dER/d6lR]/[7r7VATV2l 
given by the AdS/CFT correspondence at distances TZ — 10, 20. The AdS/CFT results are compared with the second 
order static and dynamic hydrodynamics. (7 = l/vT^-^) 



sources, we believe that this difference is due to the second order hydrodynamic coefficient, r^r 
the shear length, they are given by 



Compared to 



1.58 



V3 
c?sT 



(Boltzmann) 



0.65 



4?7/3 
c?sT 



(AdS/CFT). 



(33) 



The coefficient in unit of rj/{sT) is generically ^ 2.4 times larger in the Boltzmann equation than in 
gauge/gravity duality. As a result, the hydrodynamic solutions are more reactive to the high frequencies 
and it takes longer to reach the hydrodynamic regime, in a relative sense, in the Boltzmann equation. 
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